"Dilute" excitons in a double layer system: single-exciton and mean-field approach 
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Double layer systems where one layer has electrons and the holes are in a parallel layer a distance d 
away are expected to undergo excitonic condensation at low temperature. This excitonic condensate 
is traditionally described by a many-body wavefunction that encodes the coherence between electron 
and hole bands. Here we compare the mean-field ground state in the limit of dilute electron (hole) 
density with the ground state of a single electron-hole pair in double- layer system. As the interlayer 
distance d increases, we find that the excitonic size, characterized by the width of the momentum- 
space wavefunction, also increases. By comparing the single-exciton wavefunction with the mean- 
field analysis, we determine the d-dependence of the "diluteness" of the exciton gas in a balanced 
double-layer system with given electron (or hole) density. 



I. INTRODUCTION 

The Bose-Einstein condensation of excitons in double 
layer systems realized in semiconductor hetero junctions, 
where a macroscopic number of excitons occupy a sin- 
gle quantum state, has been a subject of intense research 
over the past decade. 1 ' 2 An exciton is a metastable bound 
state of an electron and a hole. In a balanced electron- 
hole system at low electron (hole) densities n e — nh = 
ri2D, the distance between the excitons is much larger 
than the quantum exciton size and the excitons behave 
as weakly-interacting bosons. 2,3,4 In a bulk semiconduc- 
tor or a two-dimensional system, the internal state of the 
exciton is given by Hydrogenic wavefunctions 2 that re- 
sult from the Coulomb interaction Va(i") = e 2 /er. In 
each case, the wavefunction for the electron-hole separa- 
tion decays exponentially with decay length a ex /2 where 
a ex = eh 2 /e 2 m r is the quantum size of the exciton. Here 
e ~ 10 is the dielectric constant of the semiconductor 
heterojunction, r e (r/,) represents the electron (hole) po- 
sition, m" 1 = m" 1 + 1 is the (reduced) mass of exci- 
ton, and m e (nih) is the electron (hole) band mass. Note 
that for a symmetric electron-hole system that we con- 
sider in this paper, m e = rrih = 2m r , the quantum size of 



carriers is half the exciton size, ao = eH /e m e 



c/2. 



Therefore the dimensionless distance between excitons 
for a given carrier density tl^d is given by r s /2 where 
r„, defined by 7r(r s ao) 2 = 1/«2_d is the dimensionless dis- 
tance between the carriers. Since the quantum 
of the exciton is fixed by the semiconductor properties, it 
is possible to tune the interaction between excitons from 
weak (r s 3> 1) to strong (r s ~ l). 4 Note that for a bulk 
or planar system, as opposed to a double-layer system, 
the inter-exciton interaction and the formation of exci- 
ton are both governed by the same Coulomb interaction 
V A (r). 

In double-layer systems where electrons are carriers 
in the top layer and holes are the carriers in the bot- 
tom layer, the formation of an exciton is determined by 
the attractive interlayer Coulomb interaction Ve(r) = 
—e 2 /e\/r 2 + d 2 where d is the distance between the two 



layers and r denotes the two-dimensional position vec- 
tor. Since the attractive interaction is d dependent, the 
quantum size of exciton is not necessarily a ex and, in 
fact, depends on the dimensionless ratio d/ao- The in- 
teraction between the excitons, on the other hand, is also 
dependent on intralayer Coulomb repulsion V^(r). Thus, 
in double-layer systems, the diluteness of excitons is a 
function of (d/ao,r s ). In this paper, we quantitatively 
explore this issue. In the next section we present nu- 
merical solution to the single-exciton problem. We find 
that the exciton binding energy Ef, decreases as d in- 
creases and concurrently the momentum-space ground 
state wavefunction sharpens. Thus, the quantum size of 
an exciton a ex (d) increases with d. In Sec. Ill, we com- 
plement the single-exciton results with mean-field analy- 
sis of the uniform excitonic condensate ground state for 
varying carrier density tiid- By comparing the exciton 
wavefunction obtained from the Wannier approximation 
with that in Sec. II, we obtain a quantitative criterion 
for the "diluteness" of an exciton gas. We conclude the 
paper with a remarks in Sec. IV. 



II. SINGLE-EXCITON PROBLEM 

Let us start with an electron and a hole in a double- 
layer system with d = 0. The eigenstates of this prob- 
lem are obtained by solving the equivalent problem of 
a particle with mass m r in a central attractive poten- 
tial. Due to the rotational invariance in two dimen- 
sions 5 and the existence of the conserved Runge-Lenz 
vector, 6 the energy spectrum in the limit d = is de- 
pendent only on the principle quantum number n > 1, 
E n = — AE$/{2n — l) 2 where Eq = e 2 /eao is the en- 
ergy scale associated with the problem. (Note that typ- 
ical parameters e ~ 10 and ao ~ 50A imply Eq ~ 30 
meV). The corresponding normalized ground-state wave- 
function is given by ipG^) = \/8/ {*na 2 x ) exp(— 2r/a ex ) = 
\j2j (^a 2 ,) exp(— r/a ). 5 When d ^ 0, since the electron- 
hole interaction is given by Ve(r) = —e 2 /ey/r^ + d 2 , the 
differential equation for the radial component of the ex- 
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citonic wavefunction that cannot be analytically solved. 
Instead, we use the momentum-space Schrodinger equa- 
tion 



2 1.2 



2m, 



■Va(k) 



V E (\k-k'\)^ a (k , ) = E a i; a (k) (1) 



where k is the two-dimensional wavevector, ^> a (k) is the 
momentum-space eigenfunction with eigenvalue E a , and 
Vk(q) = -V A {<l)e- qd where V A (q) = 2iTe 2 /eq is the 
Fourier transform of the intralayer Coulomb interaction 
in two dimensions. We focus on Eq.(l) projected onto the 
zero angular-momentum sector and obtain the eigenener- 
gies and eigenfunctions by discretizing the integral equa- 
tion and numerically diagonalizing the resulting matrix 7 




u n Au ~ 

—Z V(u r , 

Z7T 



(2) 



Here V(u m ,u n ) is the angular- averaged dimensionless 
electron-hole interaction. Although the Hamiltonian (2) 
leads to bound and continuum states, since the exci- 
tonic internal states are only accessible at temperatures 
T > Eo/ks ~ 300 K, in the following we only discuss the 
behavior of the ground state. 

Figure 1 shows the numerically obtained ground-state 
energy of a single exciton as a function of interlayer dis- 
tance d. At d = 0, the numerical result deviates from 
the well-known analytical answer by 10%; however, we 
have verified that this difference is solely due to dis- 
cretization errors and can be systematically suppressed. 
When d/ao *C 1 first-order perturbation theory implies 
that the change in the ground-state energy is linear, 
5E G = E G (d) - E G (0) = 4E a (d/a ). At large d the exci- 
tonic binding energy is strongly suppressed; for example, 
when d/ao — 10 it is reduced to 10% of the binding en- 
ergy at d = 0. 

In Figure 2, we show the corresponding evolution of 
the ground-state wavefunction with increasing d. At 
d = the normalized wavefunction is given by ipcik) = 
y / 87ra^/(l-|-fc 2 a 2 1 ) 3 / 2 and is reproduced by our numerical 
calculations. As d increases, we see that the momentum- 
space wavefunction sharpens and shows that the single 
exciton size a ex (d) increases with d. 

These results show that a single exciton in a double- 
layer system is increasingly weakly bound, and becomes 
larger as the interlayer distance d increases. Therefore, 
although r s is a good measure of the diluteness of elec- 
trons or holes, it is not a good measure of diluteness for 
the excitonic gas. To quantify this observation, in the 
next section, we study the evolution of a uniform exci- 
tonic condensate state as a function of (d/ao,r s ). 



III. MEAN-FIELD ANALYSIS 



FIG. 1: (Color Online) Excitonic ground-state energy Ea(d) 
as a function of interlayer distance d obtained from single- 
particle Schrodinger equation. The dotted line through d — 
shows that at small d the change in the binding energy is 
linear, 5Eq = 4Eo(d/ao), as expected from first-order pertur- 
bation theory. The binding energy is strongly suppressed at 
large d. 



7.0 

6.0 

1 5.0 

U 

£ 4-0 
3.0 



Analytical d=0.00 
Numerical d=0.00 
Numerical d=1.00 
Numerical d=2.00 



A 

°0 A 
°OA 

o. 



AO, 



ka 



FIG. 2: (Color online) Ground state wave-function xpG(k) for 
a single exciton as a function of interlayer distance d. The ana- 
lytical (cross) and numerical (open square) solutions for d = 
are consistent with each other. Their momentum-space width 
indicates that the exciton size is ao- As d increases tpG(k) 
sharpens and the size of the exciton, defined by the inverse- 
width of the momentum-space wavefunction, increases. 



distance d. The Hamiltonian for such a system is a sum 
of the kinetic energy for electrons and holes, as well as 
the intralayer Coulomb repulsion V^q) and the inter- 
layer Coulomb attraction Ve(q). We use the standard 
mean- field approximation 8,9 ' 10 to obtain the mean-field 
Hamiltonian, 



To explore the uniform excitonic condensate state, we 
start with a double-layer system with electrons in the top 
layer and holes in the bottom layer that is separated by 
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FIG. 3: (Color online) Comparison of the Wannier wavefunc- 
tions 3>(p) = Ap/Ep obtained from the mean-field solutions 
for r 3 — 7 (red circle), r a — 3 (green triangle), and r 3 — 2 (di- 
amond), with the single-exciton wavefunction Vg(p) (cross) 
for d/ao = 1. As r s increases the Wannier wavefunction ap- 
proaches the single-exciton result. 

where e k (ftJ.) is the creation operator for an electron 
(hole) with two-dimensional momentum fik in the top 
(bottom) layer, £ k is renormalized electron (hole) disper- 
sion that takes into account the interlayer capacitance 
and intralayer exchange energy, and A k is the excitonic 
order parameter associated with the coherence between 
the electron and the hole bands in the two layers. We 
consider an isotropic excitonic-condensate order param- 
eter and obtain the following self-consistent mean-field 
equations 9 ' 10 

a- = Efc + Vfc-M- / VAdk-k'IXe^ekO (4) 

A fe = y^ B (|k-k'|)(ft_ k /e k ) (5) 



n 2 D = 




where e k = h 2 k 2 /2m e = h 2 k 2 /2mi l denotes the elec- 
tron (hole) band dispersion, Vc — 27re 2 n 2 £)/e is the ca- 
pacitive energy cost, and fj, is the (electron and hole) 
chemical potential determined implicitly by Eq.(6). The 
self-consistent excitonic order parameter is determined 
by (/i- k e k ) = A k /2E k , and ±E k = ±VC k + A k denote 
dispersion of the quasiparticle bands that result from 
Hamiltonian (3). We solve Eqs.(4)-(6) iteratively for a 
given (d/ao,r s ) to obtain the self-consistent order pa- 
rameter Afc and quasiparticle energy dispersion E k . To 
explore the dilute exciton limit, we recast Eq.(5) in terms 
of $(p) = Ap/Ep : and note that for A p <C £ p Eq.(5) re- 
duces to the single-exciton Schrodinger equation in mo- 
mentum space, Eq.(l). This permits a quantitative com- 
parison between the ground-state exciton wavefunction 
4>g(p) an d the Wannier-exciton wavefunction $(p). 
Figure 3 compares the Wannier wavefunction <fr(j>) at 
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FIG. 4: (Color online) Interlayer-distance dependence of the 
critical r s value obtained using the constraints 7 = 0.90 (bot- 
tom) and 7 = 0.95 (top). r sc (d) provides a quantitative way 
to characterize the diluteness of an excitonic gas by compar- 
ing the Wannier wavefunction $(fc) with the single-exciton 
solution 4>G{k). 

d/aa = 1 for different values of r s with the single-exciton 
wavefunction ijjGip)- We see that r s increases the Wan- 
nier wavefunction approaches the single-particle result, 
as expected. Note that for small r s , the Wannier exciton 
wavefunction &(p) is peaked at finite momentum because 
the excitonic order parameter A p is maximum and £ p is 
minimum at the Fermi momentum. However, as r s in- 
creases, for any given d, the peak in A p shifts towards the 
origin and so does the maximum of the Wannier wave- 
function. 

To quantify the proximity between the Wannier and 
the single-exciton approach, we consider the overlap 
7(d/ao, between the two (real) wavefunctions 

7 (d/a ,r s ) = J $-(fc)^ G (fc). (7) 

A high overlap value 7(d/ao,r sc ) ~ 1, allows us to de- 
fine a critical value of r sc (d/ao) such that for r s > r sc 
the single-exciton result provides an excellent substitute 
for the mean- field analysis. Figure 4 shows the critical 
r sc (d) obtained using 7 = 0.90 and 7 = 0.95. We see 
that for typical values of d the critical r sc scales linearly 
with d/ao. It implies, for example, that approaching the 
dilute-limit at d/ao — 3 will require reducing the carrier 
density by a factor of 5 from the corresponding value for 
the dilute limit at d/ao = 1- 



IV. DISCUSSION 

The subject of excitonic condensation in double-layer 
systems has been extensively explored in the literature; 
the properties of a single exciton in a double-layer sys- 
tem, however, have not been. In this paper, we have 
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obtained the ground-state wavefunction and the ground- 
state energy for a single exciton as a function of interlaycr 
distance d. By comparing our results with those from a 
mean-field analysis of the uniform excitonic condensate, 
we have obtained the critical value r sc (d) that is used 
to determine when the exciton gas for a given interlayer 
distance d is "dilute" . 

Our analysis provides a quantitative picture of a single 
exciton in double-layer system with d where ana- 
lytical solution for the excitonic wave-functions is not 
possible. It shows that as d increases, due to the weak- 
ened electron-hole Coulomb interaction, the exciton size 
a ex (d) increases. 

We note that the Wannier approximation for excitonic 
wavefunction is based on a mean-field analysis that usu- 



ally over-estimates 8 the excitonic order parameter . In 
particular, in double layer systems, it is known that the 
uniform condensate becomes unstable 11 when d is larger 
than a critical layer separation d c . Thus, when fluctua- 
tions around the mean-field state are taken into account, 
the critical value of r s for a given d ~ d c will change 
substantially. Our mean-field analysis is based on a uni- 
form excitonic condensate state. Due to dipolar repulsion 
between excitons, the uniform state is unstable towards 
formation of a crystalline excitonic condensate 12 in the 
region ^/r^ < d/a < r s . Since our calculations lie out- 
side this parameter range, we have focused only on the 
uniform state; the question of a "dilute" exciton limit 
in a crystalline excitonic condensate, however, remains 
open. 
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